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Abstract 

Area-weighted Dyck-paths are a two-dimensional model for vesicles attached to a wall. 
We model the mechanical response of a vesicle to a pulling force by extending this model. 

We obtain an exact solution using two different approaches, leading to a g-deformation of 
an algebraic functional equation, and a g-deformation of a linear functional equation with a 
catalytic variable, respectively. While the non-deformed linear functional equation is solved 
by substitution of special values of the catalytic variable (the so-called "kernel method"), 
the g-deformed case is solved by iterative substitution of the catalytic variable. 

Our model shows a non-trivial phase transition when a pulling force is applied. As soon 
as the area is weighted with non-unity weight, this transition vanishes. 

1 Introduction 

With the growing experimental ability to probe the behaviour of a single polymer under the 
influence of a stretching force \13\ 112] , there has been an effort by theoretical work to explore 
models that cover the possible scenarios encountered. For example, various directed path models 
of a polymer with one end, or both ends, grafted to a surface and being pulled at the end or 
in the middle, with the addition of an attractive contact interaction with the surface, have 
been studied in two and three dimensions |22[ [2"51 [TJ [T5]. Numerical work on self-avoiding 
walks models has also confirmed and expanded upon the predictions of these directed models 
|16} \T7\ [20] I21j . In the continuum limit, this problem is related to the area under a Brownian 
curve [19]. 

Another type of model, that of a lattice polygon that has been weighted by its length and 
area, can rather be used to model biological vesicles JTHJ, |2] [91 [231 OH]. For example, the effect 



1 



of a restricted geometry such as a slit has recently been analysed [25]. Here, it is also natural 
to consider the effect of a pulling force of the surface of the vesicle. 

Pulling membrane- anchored microspheres using optical traps is a technique commonly used 
in experimental settings [28]. When preparing the sample, microspheres can bind anywhere on 
the membrane, and due to small anisotropies in the experimental setup perfect vertical pulling 
is difficult to achieve. 

Hence, we study a simple directed model of a two-dimensional vesicle that is being pulled 
away from a wall. The anisotropy is modelled by a pulling force with a small horizontal com- 
ponent, which leads to the force acting on the right-most peak of the vesicle. Only the vertical 
component of the force contributes to the configurational energy of the vesicle. 

The model is of interest physically as well as mathematically. It illustrates the solution 
of functional equations with catalytic variables under a (/-deformation, when the solution gets 
expressed in terms of g-series rather than in terms of algebraic functions. 



Figure 1: An area- weighted Dyck path with 18 steps and area of size 11, which is being pulled 
away from the surface at its right-most peak with force /. The height h of the right-most peak 
in this example is 4. 

Dyck paths are directed walks on Z 2 starting at (0, 0) and ending on the line y = 0, which 
have no vertices with negative y-coordinates, and which have steps in the (1, 1) (up-step) and 
(1, —1) (down-step) directions [29]. Given a Dyck path tt, we define the length n(n) to be half 
the number of its steps, and the area m(ir) to be the sum of the starting heights of all up-steps 
in the path. This is sometimes also called the rank function of a Dyck path [8j |27], and is 
equivalent to the number of diamond plaquettes under the Dyck path. We further define the 
height /i(vr) of the right-most peak as the length of the final descent of the Dyck path, i.e. the 
maximal number of the last consecutive down-steps. An example of a Dyck path is given in 



2 The Model 




Fig. EQ 
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Let Q be the set of Dyck paths, and define the generating function 

G(t, q,\) = Y J t n ^q m ^\ hM . (2.1) 
Tree 

The generating function for area- weighted Dyck paths |10| [7] has previously been studied. Dyck 
paths weighted by height of the right-most peak are related to Ballot paths, which have been 
previously studied as a model of pulled polymers [22] . Here we extend these works by combining 
and extending these models. 

The key to obtaining a solution of this model is a combinatorial construction of Dyck paths 
which leads to a functional equation for the generating function G. We shall present two different 
constructions leading to two structurally different functional equations. 

3 g-deformed Algebraic Functional Equation 

To obtain the first functional equation, we decompose Dyck paths by considering the right-most 
up-step starting at height zero. This leads to a decomposition of a Dyck path in terms of two 
smaller Dyck paths. More precisely, except for the zero-step path, this decomposes a Dyck path 
uniquely into a Dyck path, followed by another Dyck path which is bracketed by up and down 
steps, as indicated in Figure [2j 



Figure 2: A schematic decomposition leading to Eqn. (3.1): A non-zero Dyck path decomposes 
uniquely into a Dyck path, followed by another Dyck path which is bracketed by up and down 
steps. 

This naturally leads to a functional equation which is quadratic in the generating functions, 

G(t, q, A) = 1 + XtG{t, q, l)G(qt, q, A) . (3.1) 

Note that in the product XtG(t,q, l)G(qt, q, A) the left Dyck path no longer contributes to the 
height of the right-most peak, leading to the factor G(t,q, 1), and the bracketing of the right 
Dyck path by up and down steps increases the area proportional to the length of the path, as 
well as increases height and length by one, leading to the factor XtG(qt, q, A). 

3.1 The Case q = 1 



The functional equation (3.1) reads now 

G(t, 1, A) = 1 + XtG(t, 1, l)G(t, 1, A) , (3.2) 
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which after substitution of A = 1 can be readily solved to give 

1 - y/T=M 2 



G(t,M) 



(3.3) 



2t 1 + ^T^U 

which is the well-known generating function for Catalan numbers. Substitution of this result 
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back into Eqn. (3.2) gives 



G(t,l,X) 



For comparison with the result in Eqn. (3.11), note that a series expansion in A gives 



G(t,l,X)=Y / \ l t l G(t,l,l) 1 



(3.4) 



(3.5) 



3.2 The case q ^ 1 



Substituting A = 1 into the functional equation (3.1) gives 



The Ansatz 



G(t,q,l) = l + tG(t,q,l)G(qt,q,l) 
H(qt,q) 



G(t,q,l) 



H(t,q) 



gives a linear functional equation for H(t, q), 

H{qt,q)=H(t,q)+tH{q 2 t, 



(3.6) 
(3.7) 

(3.8) 



Inserting H(t,q) = Y^=o c n{l)t n into (3_8) gives a recurrence q n c n (q) = c n (q) + q 2{ ~ n ^c n -i(q) 
for the coefficients, which can be readily solved by iteration. Specialising to co(q) = 1, we find 



~ (-i) k a k (k-i) 



(3.9) 



where we have used the q-product notation 



k-i 



(t;q) h = H(l-tqi). 

0=0 



(3.10) 



Here H(t, q) = o(j)i(—; 0; q, —t) is a basic hypergeometric series [TT] also referred to as the g-Airy 
function |14j . 

Once G(t,q, 1) is known and substituted back into Eqn. (3.1), we can solve for G(t,q,X) by 
iteration and obtain 

G(t.,.A) = f A¥,(«^>. (3.11) 



1=0 



This iteration makes sense in a ring of formal power series. Alternatively, one can argue that 
the iteration converges for \q\ < 1. 

We note that this solution also admits a direct combinatorial interpretation. The Z-th term in 



the sum of (3.11 ) corresponds to the concatenation of / Dyck paths separated by single up-steps, 
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Figure 3: The term I = 3 in the sum of (3.11) corresponds to the concatenation of three Dyck 
paths separated by single up-steps, followed by a hook formed of one up-step and three down- 
steps. 



followed by a final hook formed by one up-step and I down-steps. The generating function for a 
Dyck path starting at height I is given by G(q l t, q, 1), and the final descent of length I increases 
the by ( 2 ) , leading to an overall weight of 



G(t, q, l)tG(qt, q,l)t-.... G(q l t, q, l)t\ l q® = X l t l q^) ^jr4 . 

Figure [3] indicates this for the case 1 = 3. 

Finally, we note that the limit q — > 1 can be obtained as follows. As 

G(t, 1,1) = lim G(t, q, 1) = lim > 

H(t,q) 



we deduce that for fixed value of k, 



G(t, 1,1) = lim 



H(q k+1 t,q) 



and thus 



q^l H(qH,q) 

v H(q l t,q) , 
q^i H(q,t) 



(3.12) 



(3.13) 

(3.14) 
(3.15) 



Therefore Eqn. (3.11) reduces to the result in Eqn. (3.5) as q — > 1. 



4 (/-deformed Linear Functional Equation 

To obtain the second functional equation, we consider how Dyck path configurations get changed 
by the insertion of a small pattern. 

The right-most descent, together with the preceding up-step forms a hook. Inserting an up- 
step followed immediately by a down-step in the right-most descent creates a new Dyck path 
with a new hook. If the descent has length h, there are h + 1 possibilities of insertion. 

The consequence of this insertion is the addition of a column to the Dyck path as indicated 
on the left-hand-side of Figure |4} 

The functional equation corresponding to this construction is given by 

G(t, q,X) = l+ j^Git, q, 1) - ^G(i, q, qX) . (4.1) 
1 — qX 1 — qX 
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Figure 4: A schematic construction leading to Eqn. (4.1): Inserting an up-step followed im- 
mediately by a down-step in the right-most descent creates a new Dyck path with a new hook 
(left). When adding hooks of all possible lengths, one needs to correct for overhangs produced 
by hooks that are too long (right). 

The term i^t\ G(t, q, 1) corresponds to the addition of a column of arbitrary height to 
the Dyck path. The ensuing overcounting is compensated for by subtraction of the term 
l-q\ G(t, q, qX), which corresponds to the removal of the columns that produce invalid con- 
figurations, as indicated on the right-hand-side of Figure [4} 

Alternatively, the terms on the right-hand side of {14. 1 h can be understood by considering 



G(t,q,X) = Y / X h G h (t,q), (4.2) 



h=0 

where Gh(t,q) are partial generating functions for fixed height h. The only Dyck path with a 
hook of height h = is a zero-step path. A Dyck path with a hook of height h > 1 can only have 
been constructed from a Dyck path with a hook of height at least h — 1. By the construction, 
its area increases by an amount h — 1, and we thus find that for h > 1 

oo 

G h (t,q)=q h ~ 1 Gk(t,q). (4.3) 

k=h-l 



Summing over the partial generating functions leads to the functional equation (4.1). 



4.1 The Case q = 1 



The functional equation (4.1) now reads 



G(t, 1, A) = 1 + ^G(t, 1, 1) - ^jG(t, 1, A) . (4.4) 
Here, substitution of A = 1 is no longer useful, and one solves the functional equation using the 



kernel method [26] . Writing Eqn. (4.1) as 



K(t, X)G(t, 1, A) = 1 + Y^G(t, 1, 1) , (4.5) 



where the kernel K(t, A) is given as 



K(t,X) = l + ^, (4.6) 
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we observe that the left-hand side of Eqn. (4.5) vanishes if A = Ao is chosen such that K(t, Xq 
vanishes (provided G(t, l,Ao) does not diverge). This then implies that 

Ao-1 



Ant 



(4.7) 



and substituting this result back into Eqn. (J45J) leads to an expression for G(t,l,X). The 
appropriate root of the kernel is given by 

1 - VT^tt 



Ao 



2t 



(4.8) 



Using this value, we recover the expressions (3.3) for G(t, 1,1) and (3.4) for G(t, 1, A) given 
above. 



4.2 The Case q ^ 1 



If q 1, we can solve the functional equation (4.1) by iterative substitution. Substituting Xq 



for A in (4.1), we find 



G(t,q,q k X) = l + 



q k Xt 



We can therefore express G(t,q,q X) in terms of G(t,q,q + X) and use this to iterate. As a 



2k+l \2 f 
k+1 



(4.9) 



result, we find 



Defining 



this simplifies to 



G(t,q,X) = J2 



{-X 2 t) k q k 
(<?A;g); 



fc=0 



I(t,q,y) = J2 



k=0 



1 



n k\f 

1 + d^TTA G( ^' r 

(_ t J*g*(fc"l) 

(y;?)fc 



G(t, g, A) = /(gA 2 i, g, gA) - - [/(A 2 t, g, qX) - l] G(t, g, 1) 



(4.10) 



(4.11) 



(4.12) 



Note that I(t,q,y) = i</>2(g; 0, y; g, — t) is also expressible as a basic hypergeometric function. 
Moreover, it a deformation of H(t, g), in that I(t, g, g) = H(t, g). 



Substituting A = 1 in Eqn. (4.12) gives now immediately 

H(qt,q) 



G(t,q,l) 



(4.13) 



a result already obtained from the g-algebraic functional equation (3.6) above. Substituting this 



result back into Eqn. (4.12) leads to the main result of this section, namely 



G(t, q, A) = I{qX\ q, qX) - ± [/(A 2 t, q, qX) - l] . 



(4.14) 



Note the rather different structure of the equivalent solutions (3.11) and (4.14). The former is 



in terms of an infinite sum, which however readily allows for extraction of the partial generating 
function for fixed height h. The latter is more compact as it is written in terms of a finite 
number of special functions. We note that the equivalence of both solutions implies the identity 

ns __, , , _, . n . . __. . 1 



X l t l q^H(q l t, q) = I(qX 2 t, q, qX)H(t, q) - j [l(X 2 t, q, qX) - l] H(qt, q) . 

1=0 



(4.15) 
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5 Analysis of behaviour 



0.3 



0.2 



0.1 



3 
X 



Figure 5: A plot of the radius of convergence t c (X) = i c (A, 1) in the polymer length variable t 
of the generating function when q = 1 as a function of the force variable A. Note that there is a 
change of behaviour at A = 2 with constant t c (\) = 1/4 for A < 2. 

Let us begin by considering the case of q = 1, that is no extra area weighting or, equivalently, 
osmotic pressure. The first important point to make is that previous models of pulling on the 
end of a polymer without any other interaction with a surface give the result that as soon as a 
force is applied the end of the polymer changes its scaling behaviour from random-walk scaling 



(h) ~ n 1 ' 2 



for / = (A = 1) to ballistic scaling 



(h) 



n 



(5.1) 



(5.2) 



for / > (A > 1). This is reflected in the fact that the radius of convergence of the associated 
generating function is constant for A < 1 and decreases continuously for A > 1. As a model of 
pulling on a vesicle attached to a surface, our model has the feature that the pulling occurs on 
the final "peak" , the location of which along the surface can change with the pulling itself. As 
such, there is a non-trivial critical point as a function of A in our model, with a critical value 

A = A c = 2 . (5.3) 

A plot of the radius of convergence of the generating function (3.4) is shown in Fig. [5j 
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Note that the reduced free energy k(X, q) of the surface is simply related to the radius of 
convergence t c (X, q) in the usual way via 



= - log t c (\,q) . 



(5.4) 



For q < 1 the generating function (3.11 ) can be seen to be convergent except where H(t, q) 



0. Crucially this is independent of A! This implies that for q < 1 there is no phase transition 
as a function of A, since n(X,q) is independent of A for q < 1. Physically this means that the 
osmotic pressure always dominates the pulling force for q < 1: (h) will be bounded for all A. 

Interestingly, if one considers the limit q — > 1 for A > 2 one is left with the problem that the 
free energy "jumps" at q = 1 as indicated in Figure [6j The resolution of this apparent problem 
is that for q > 1 one needs to consider the free energy per unit area rather than the free energy 
per unit length of the surface. The limiting free energy per unit area is a continuous function 
of q and A for all A, being equal to zero for q < 1. 




Figure 6: A plot of the radius of convergence t c (X,q) in the polymer length variable t of the 
generating function for general q as a function of the force variable A. Note that t c (X,q) is 
independent of q for q < 1 and hence that there is a jump of i c (A, q) as q approaches one from 
below. 



6 Conclusion 

We have provided the exact solution of a simple model of a vesicle attached to a surface subject 
to a pulling force on the vesicle. We have used two methods of solution which illuminate the 



9 



generalisation of the well-known "kernel" method for the solution of linear functional equations 
using a "catalytic" variable and result in a non-trivial identity. The solution is expressed in 
terms of well-studied ^-series and can be analysed simply. Our model has a non-trivial phase 
transition when no osmotic pressure is applied. 
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